Strong frustration due to competing ferromagnetic and antiferromagnetic 
interactions: magnetic properties of M(VO)2(P04)2 (M = Ca and Sr) 
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We present a detailed investigation of the magnetic properties of complex vanadium phosphates 
M(VO)2(P04)2 (M — Ca and Sr) by means of magnetization, specific heat, "^^P NMR measure- 
ments, and band structure calculations. Experimental data evidence the presence of ferromagnetic 
and antiferromagnetic interactions in M(VO)2(P04)2, resulting in a nearly vanishing Curie-Weiss 
temperature &cw < 1 K that contrasts with the maximum of magnetic susceptibility at 3 K. Specific 
heat and NMR measurements also reveal weak exchange couplings with the thermodynamic energy 
scale Jc = 10 — 15 K. Additionally, the reduced maximum of the magnetic specific heat indicates 
strong frustration of the spin system. Band structure calculations show that the spin systems of the 
M(VO)2(P04)2 compounds are essentially three-dimensional with the frustration caused by com- 
peting ferromagnetic and antiferromagnetic interactions. Both calcium and strontium compounds 
undergo antiferromagnetic long-range ordering at Tm = 1-5 K and 1.9 K, respectively. The spin 
model reveals an unusual example of controllable frustration in three-dimensional magnetic systems. 
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PACS numbers: 75.50.Ee, 75.40.Cx, 71.20.Ps, 75.30.Et 



I. INTRODUCTION 

Low-dimensional and/or frustrated spin systems are a 
subject of thorough studies due to their variety of un- 
usual ground states and quantum phenomena. ^'^ One of 
the most striking effects in these systems is the forma- 
tion of a spin-liquid, a dynamically disordered ground 
state with short-range magnetic correlations. Spin liq- 
uids are caused by quantum fluctuations that impede or 
even suppress long-range magnetic ordering 

Quantum fluctuations are enhanced in systems with 
low spin value, low dimensionality, and/or magnetic frus- 
tration. At present, numerous spin models (and the re- 
spective structural types) for the frustrated spin systems 
are known. '^''^ Most of these systems are geometrically 
frustrated magnets, since the frustration is caused by the 
topology of the lattice. However, several models reveal 
frustration due to a specific topology of magnetic interac- 
tions, while the lattice itself does not prevent the system 
from long-range spin ordering. An important advantage 
of the latter systems is the possibility to vary the degree 
of frustration by tuning individual exchange couplings. 
Hence, the formation of different ground states within 
one model and the access to quantum critical points 
are possible. Below, we focus on the respective models 
and briefly review both one-dimensional (ID) and two- 
dimensional (2D) cases. 

The simplest ID lattice is a uniform chain. To frus- 
trate such a system, one should introduce next-nearest- 
neighbor (NNN) interaction J2. Typically, J2 is antiferro- 
magnetic and competes with either ferromagnetic (FM) 
or antiferromagnetic (AFM) nearest-neighbor (NN) in- 
teractions Ji. The resulting physics strongly depends 
on the sign of Ji. Thus, for Ji > (i.e., for AFM NN 
coupling) the system has a dimerized ground state for 



J2/J1 > 0.241 that may be further stabilized by lat- 
tice distortion as experimentally observed in the spin- 
Peierls compound CuGeOai^i^ If Ji < 0, hehcal mag- 
netic ordering is formed for J2/J1 < —0.25, and the 
pitch angle of the spiral depends on the J2/J1 ratio. 
Recently, a number of frustrated chain systems with 
J2/J1 < have been studied,'' and even the vicinity of 
the quantum critical point at J2/J1 = —0.25 was ac- 
cessed experimentally)^ In these systems, the possibility 
to tune the pitch angle by varying the exchange cou- 
plings may be particularly attractive due to another in- 
triguing property - the magnetoelectric effect observed 
in LiCusOs and UCuV04^'^^^^^^^ 

In two dimensions, we will consider the square lattice. 
Frustration can be introduced by diagonal NNN inter- 
actions (J2) that compete with NN interactions Ji run- 
ning along the side of the square. The resulting model 
is known as a frustrated square lattice (FSL) or Ji — J2 
model. The phase diagram of the FSL has been exten- 
sively studied theoretically (see, e.g., Refs. [l3lfl3ll5lll6l ). 
It shows three ordered phases (ferromagnet, Neel anti- 
ferromagnet, and columnar antiferromagnet) as well as 
two critical regions near the quantum critical points at 
J2/J1 — ±0.5. Basically, these regions are supposed to 
reveal spin-liquid ground states, although there is also 
a proposition of a nematic order at J2/J1 = — O.Sfii Yet 
both the conjectures are not verified experimentally, since 
most of the FSL systems studied so far show columnar 
magnetic ordering and do not fall to the critical regions. 
Nevertheless, the detailed investigation of the respective 
compounds was highly helpful and disclosed important 
methodological aspects of studying the frustrated spin 
systems and understanding their physics J^ii^i^Sii^ 

Based on the above representative examples, one may 
expect that similar changes of the ground state due to the 
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alteration of the exchange couphngs could be realized in 
other frustrated models - either in one, two, or even three 
dimensions. This consideration and the recent investiga- 
tion of an interesting FSL system, the complex vana- 
dium phosphate Pb2VO(P04)2 (Ref. [12), motivated us 
to start a systematic research of the related compounds. 
Below, we present the study of M(VO)2(P04)2 (M = Ca, 
Sr) - two representatives of the same compound family 
that reveal three-dimensional (3D) spin systems with the 
frustration caused by the specific topology and magni- 
tudes of the exchange interactions. In the following, we 
show clear evidence of the frustration in the results of 
several experimental techniques: magnetization, specific 
heat, and nuclear magnetic resonance (NMR) measure- 
ments. In order to understand the microscopic origin of 
the frustration, we use band structure calculations. 

The outline of the paper is as follows. In Sec. |TT1 we 
analyze the crystal structure of the M(VO)2(P04)2 com- 
pounds focusing on possible paths of the exchange inter- 
actions. Methodological aspects are reviewed in Sec. IIIII 
Section llVI presents experimental results including mag- 
netization, specific heat, and '^^P NMR. In Sec. El we 
report the results of band structure calculations and con- 
struct a microscopic model of the exchange interactions 
in Ca(VO)2(P04)2. Experimental and theoretical results 
are compared and discussed in Sec. IVII followed by our 
conclusions. 



II. CRYSTAL STRUCTURE 

Complex vanadium phosphates Ca(VO)2(P04)2 
(Ref.lH) and Sr(VO)2(P04)2 (Ref. |2j) are isostructural 
and crystallize in a face-centered orthorhombic unit cell 
(space group Fdd2, Z = 8). The structure presents 
one vanadium site that shows distorted octahedral 
coordination. The V+^Og octahedra share corners to 
form chains running along the [101] and [101] directions. 
The chains are joined into a 3D-framework by PO4 
tetrahedra (Fig. [1]). 

Magnetic interactions in insulating transition-metal 
compounds mainly depend on two factors: the orbital 
state of the magnetic cation and the possible paths for 
superexchange interactions. In the following, we will dis- 
cuss both of these factors. 

The orbital state of the magnetic cation is determined 
by the local environment. In the M(VO)2(P04)2 com- 
pounds, vanadium has distorted octahedral coordination 
with one short bond (see bottom part of Fig. [T]). This 
type of environment is typical for V^** (electronic config- 
uration 3d^) and gives rise to a non-degenerate ground 
state with half-filled dxy orbital (see Refs. [2^ andf2^ and 
Sec. |Vl. The dxy orbital is located in the plane perpen- 
dicular to the short V-0 bond. Short bonds are aligned 
along the structural chains, therefore the overlap of the 
half-filled orbital with p orbitals of the two bridging oxy- 
gen atoms is very weak, and V-O-V is not an efficient 
path for superexchange interactions. On the other hand, 




FIG. 1: (Color online) Projection of the Ca(VO)2(P04)2 crys- 
tal structure along the [101] direction (upper panel) and the 
chain of corner-sharing VOe octahedra (bottom panel). In 
the upper panel, the arrows denote different paths of superex- 
change interactions running via PO4 tetrahedra. To simplify 
the figure, calcium atoms are not shown. In the bottom panel, 
thick (yellow) lines indicate short V-0 bonds that are aligned 
along the direction of the chain. 



the dxy orbital may overlap with the p orbitals of four 
other oxygen atoms providing more complicated V-0- 
P-O-V paths for the superexchange. 

The chains of the VOg octahedra in the structures of 
the M(VO)2(P04)2 compounds are not parallel, and the 
V-P-0 framework looks quite complicated. In order to 
simplify the situation, one should take into account two 
considerations: (i) Every vanadium atom is connected 
to four different PO4 tetrahedra. (ii) Every PO4 tetra- 
hedron is connected to four different vanadium atoms. 
Thus, every vanadium atom may have 3x4= 12 dif- 
ferent V-O-P-O-V connections at the most. In fact, 
we find only ten of them, since two connections (namely, 
that corresponding to J3) are provided by two tetrahedra 
simultaneously (see Fig.[T]). Some of the connections are 
equivalent by symmetry and hence give rise to identical 
exchange interactions. Finally, we find that every vana- 
dium atom can take part in five different superexchange 
interactions running via PO4 tetrahedra, and there are 
two interactions of each type per vanadium atom. We 
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FIG. 2: (Color online) Reciprocal magnetic susceptibility (1/x) vs. temperature (T) measured at /loH = 0.5 T down to 0.4 K 
for Ca(VO)2(P04)2 (left panel) and Sr(VO)2(P04)2 (right panel). The solid lines are best fits of the data to Eq. In the 
insets, x vs. T is plotted in the logarithmic scale to focus on the T/^ax a^nd Tjv. 



denote these interactions and corresponding hoppings by 
Ji — J5 and ti — is, respectively (see Fig. [1]). Ji runs 
along structural chains, J2 runs between parallel struc- 
tural chains, while J3 — J5 provide couplings between 
non-parallel chains. 

According to a recent study of another vanadium phos- 
phate, Ag2VOP207 (Ref. i25|), exchange couphngs via 
PO4 tetrahedra can hardly be estimated using struc- 
tural data only. Such couplings depend on numerous 
geometrical parameters, and slight structural alterations 
may give rise to a huge change of the exchange inte- 
grals. Therefore, we do not attempt to construct a spin 
model based on qualitative and empirical considerations. 
Instead, we proceed to the results of band structure 
calculations and estimate individual exchange couplings 
(Sec. |V|. One should be aware that the spin systems of 
the M(VO)2(P04)2 compounds are intricate and quite 
complicated even for the computational study due to the 
weakness of the exchange couplings and the presence of 
numerous NN and NNN exchange pathways. Therefore, 
the computational estimates are less accurate as com- 
pared to that in other spin-1/2 systems (see, e.g., Refs. 
H, [13, and I25I). Nevertheless, we succeed to reveal the 
leading exchange couplings and suggest a plausible sce- 
nario of the magnetic frustration in the M(VO)2(P04)2 
compounds (see Sec. IVI|) . 



III. METHODS 

Polycrystalline samples of M(VO)2(P04)2 were pre- 
pared by solid state reaction technique using CaCOs 
(99.9%), SrCOa (99.999%), VO2 (99.99%), and 
(NH4)2HP04 (99.9%) as starting materials. The syn- 
thesis involved two steps. In the first step, the inter- 
mediate compounds MP2O6 were prepared by firing the 



stoichiometric mixtures of MCO3 and (NH4)2HP04 at 
850 °C for 1 day in air with one intermediate grind- 
ing. In the second step, the intermediate products 
were mixed with VO2 in the appropriate molar ratio 
and heated for 4 days in dynamic vacuum (10""" mbar) 
with several intermediate grindings and pelletizations. 
The annealing temperatures were 900 and 850 °C for 
Ca(VO)2(P04)2 and Sr(VO)2(P04)2, respectively. The 
resulting samples were single-phase as confirmed by X- 
ray diffraction (STOE powder diffractometer, CuKq, ra- 
diation). Lattice parameters were calculated using a 
least-squares fit procedure. The obtained lattice pa- 
rameters for Ca(VO)2(P04)2 [a = 11.774(2) A, 6 = 
15.777(3) A, and c = 7.163(1) A] and Sr(VO)2(P04)2 
[a = 12.008(2) A,b= 15.924(3) A, and c = 7.207(1) A] 
are in agreement with the previously reported valuesi^!^ 

Magnetization data were measured as a function of 
temperature and field using a SQUID magnetometer 
(Quantum Design MPMS). Heat capacity data were col- 
lected with Quantum Design PPMS on pressed pellets 
using the relaxation technique. All the measurements 
were carried out over a large temperature range (0.4 K 
< T < 400 K) in magnetic fields fioH up to 14 T. The 
low-tcmperature measurements were done partly using 
an additional "^He setup. 

The NMR measurements were carried out using pulsed 
NMR techniques on ^^P (nuclear spin 1=1/2 and gy- 
romagnetic ratio 7/27r = 17.237 MHz/Tesla) nuclei in a 
large temperature range (1 K < T < 300 K). We have 
done the measurements at a radio frequency of 70 MHz 
which corresponds to an applied field of about 4.06 T. 
Low-temperature NMR measurements were partly done 
using a '^He/^He dilution refrigerator (Oxford Instru- 
ments) with the resonant circuit inside the mixing cham- 
ber. Spectra were obtained either by Fourier trans- 
form (FT) of the NMR echo signals or by sweeping the 
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field at fixed frequency of 70 MHz. The NMR shift 
K{T) — [ly (T) — z/ref] /vrei was determined by measuring 
the resonance frequency of the sample {i^ (T)) with re- 
spect to a standard H3PO4 solution (resonance frequency 
fref)- The '^^P spin-lattice relaxation rate (l/Ti) was 
measured by using either a 180° pulse (inversion recov- 
ery) or a comb of saturation pulses. 

Scalar-relativistic band structure calculations were 
performed using the full-potential local-orbital scheme 
(FPL07. 00-27)^ and the parametrization of Perdew and 
Wang for the exchange and correlation potential.^'* A k- 
mesh of 1728 points within the first Brillouin zone (510 in 
the irreducible part for the space group Fdd2) was used. 
Convergence with respect to the fc-mesh was carefully 
checked. 

First, a local density approximation (LDA) calculation 
was done using the full symmetry of the crystal struc- 
ture (space group Fdd2). The results of this calculation 
enabled the choice of the orbital states relevant for the 
magnetic interactions. The respective bands were ana- 
lyzed using a tight-binding (TB) model, and the resulting 
hoppings were applied to estimate the antiferromagnetic 
contributions to the exchange integrals. Such an analysis 
accounts for all reasonable paths of superexchange and 
provides a reliable basis for the microscopic model of the 
magnetic interactions. However, LDA does not give any 
information about ferromagnetic couplings in the system 
under investigation. 

To account for ferromagnetic interactions in the mi- 
croscopic model, we performed LSDA-I-C/ calculations 
and compared total energies for four different patterns of 
spin ordering. Exchange couplings in the M(VO)2(P04)2 
compounds are very weak; therefore one has to use the 
smallest unit cell (i.e., the primitive cell - one- fourth 
of the face-centered unit cell) in order to achieve the 
highest possible accuracy. The primitive cell is triclinic 
(space group PI) and includes four inequivalent vana- 




dium atoms, hence a number of different spin configura- 
tions can be formed. 

Following the results oin Ref. [H, we employed sev- 
eral values of Ud (Coulomb repulsion parameter of the 
LSDA+C/ method) in our calculations. According to the 
previous study^^ J7d = 6 eV provides a reasonable de- 
scription of the octahedrally coordinated V"^* within the 
FPLO calculations. Yet one should keep in mind that the 
optimal value of Ud depends on numerous factors (basis 
set,^^ local environment of the transition-metal cation, 
and objective parameters), and the unique choice of Ud 
remains a subtle issue. Below, we use several represen- 
tative values of Ud (4, 5, and 6 eV) and carefully check 
their effect on the results. The exchange parameter of 
LSDA-f [/ was fixed at J = 1 eV since usually it has mi- 
nor importance compared to Ud, especially for the 2>d^ 
configuration. 



IV. EXPERIMENTAL RESULTS 

A. Bulk susceptibility and magnetization 

Bulk magnetic susceptibilities (x) of the 
M(VO)2(P04)2 compounds are presented in Fig. [2l 
With decreasing temperature, the susceptibility in- 
creases in a Curie- Weiss manner and passes through a 
broad maximum at TjJax — 3 K. The maximum is char- 
acteristic of low-dimensional spin systems and indicates a 
crossover to a state with antiferromagnetic correlations. 
The change in slope observed at T/v = 1.5 K and 1.9 K 
[for Ca(VO)2(P04)2 and Sr(VO)2(P04)2, respectively] 
can be attributed to the onset of long-range magnetic 
ordering. 

To fit the bulk susceptibility data at high tempera- 
tures, we use the expression: 



X = Xq + 



c 



T + 6'cw 



(1) 



FIG. 3: Isotliermal magnetization of tlie M(VO)2(P04)2 com- 
pounds measured at 1.8 K. 



where xo is the temperature-independent contribution 
that accounts for core diamagnetism and Van Vleck para- 
magnetism, while the second term is the Curie- Weiss law 
with the Curie constant C — NAfJ^'^ff/Sks- The data 
above 20 K were fitted with the parameters listed in Ta- 
ble [D The resulting effective moments /i^^j are in rea- 
sonable agreement with the spin-only value of 1.73 /is, 
while the Curie- Weiss temperatures for both the 

compounds are rather small as compared to T^^x 
use the superscript x to denote the values related to the 
analysis of the bulk susceptibility data) . 

Basically, T^ax characterizes the energy scale of the 
magnetic interactions, and 6cw is a linear combination 
of all the exchange integrals.'^- Therefore, the reduction 
in ^cw as compared to TjJax implies the presence of both 
FM and AFM interactions in the system under investiga- 
tion. Moreover, Tjv is considerably lower than T^ax' ^^'^ 
this effect may be caused by either low dimensionality 
and/or frustration of the spin system. As we will show 
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TABLE I: The parameters of the Curie- Weiss fit [Eq. ^] of the bulk susceptibility and the NMR shift data (marked by the 
superscripts x a-nd K, respectively), xo is the temperature-independent contribution, fics is the effective magnetic moment, 
and 6cw is the Curie- Weiss temperature. The listed standard deviations originate from the least-squares fitting. 



Sample 


Xo (10 ^ emu/mole) 




egw (K) 




^cw (K) 


Ca(VO)2(P04)2 
Sr(VO)2(P04)2 


2.1(5) 
-1.5(5) 


1.630(1) 
1.653(1) 


-0.4(3) 
1.0(1) 


1.71(6) 
1.70(5) 


-0.3(1) 
-0.4(1) 



below (Sec. E]), the spin systems of the M(VO)2(P04)2 
compounds are three-dimensional. Thus, the reduction 
in Tn as compared to T^^^ indicates magnetic frustration 
in M(VO)2(P04)2. 

Field-dependent magnetization data for the 
M(VO)2(P04)2 compounds are shown in Fig. [3l At low 
fields, the curves reveal linear behavior, while a positive 
curvature is observed at higher fields. Further increase 
in the field results in the saturation at noHg ~ 8 T 
and 11.5 T for Ca(VO)2(P04)2 and Sr(VO)2(P04)2, 
respectively. The saturation magnetization (Ms) is 
about 0.9 /is/mole, i.e., slightly below the expected 
value of 1 ^B- The underestimate of may be 
caused by non-magnetic impurities in the samples under 
investigation. This explanation is further supported by 
the slight reduction in and the magnetic entropy 
(see Sec. UVBl) . 

The saturation field shows the energy difference be- 
tween the ground state and the fully polarized (ferro- 
magnetic) state of the system. Thus, the values of Hs 
can be used to estimate exchange couplings, if the ground 
state of the system and the leading exchange interactions 
are known. We will further discuss this point in Sec. IVII 
and employ the Hs values to get quantitative information 
about exchange couplings in the systems under investi- 
gation. 



B. Specific heat 

The specific heat (Cp) results at zero field are shown 
in Fig. 131 At high temperatures, Cp is completely dom- 
inated by the contribution of phonon excitations; there- 
fore, both the compounds reveal similar Cp{T) curves. 
With decreasing temperature, Cp starts to increase be- 
low 12 — 15 K indicating that the magnetic part of the 
specific heat (Cmag) becomes prominent. With a further 
decrease in temperature, Cp (T) shows a broad maximum 
at around 3 K and 4 K due to the correlated spin exci- 
tations and a sharp peak at Tjv = 1.5 K and 1.9 K as- 
sociated with the long-range magnetic ordering [the val- 
ues are given for Ca(VO)2(P04)2 and Sr(VO)2(P04)2, 
respectively] 

In order to get a quantitative estimate of Cmag, the 
phonon part Cphon was subtracted from the total mea- 
sured specific heat Cp. The general pro cedure is sim- 
ilar to that reported in Refs. |32 and l33l The phonon 
part was estimated by fitting Cp{T) at high temperatures 



(15 K < T < 200 K) with a sum of Debye contributions. 
The additional term A/T'^ accounted for the magnetic 
contribution, and the final fit was performed using the 
equation: 

„ (2) 

where R = 8.314 J/molc K is the gas constant, q'^^ are 
the characteristic Debye temperatures, and Cn are inte- 
ger coefficients indicating the contributions of different 
atoms (or groups of atoms) to the specific heat. The 
phonon contribution was extrapolated down to 0.4 K 
and subtracted from the measured Cp{T). The relia- 
bility of the whole procedure was justified by integrating 
Cmag/T. The resulting magnetic entropies are 5.30 and 
5.58 J/mole K [for Ca(VO)2(P04)2 and Sr(VO)2(P04)2, 
respectively] consistent with the expected value of i? In 2. 

High-temperature magnetic contribution A/T"^ is the 
lowest-order term in the high-temperature series expan- 
sion for the specific heat. According to Ref. [30l 

i 

where integers Zi indicate the number of interactions 
Ji for a single magnetic atom (i.e., coordination num- 
ber for the interactions of type Ji). The structures of 
M(VO)2(P04)2 yield = 2 for any i (see Sec. Ill, and 
we find A = {3R/16)J^ with the thermodynamic energy 
scale of the exchange couplings defined as Jc = ^/J2i Jf ■ 
Using the experimental values A = 142 and 147 J 
K/mole, we estimate = 9.6 K and 9.8 K for the 
calcium and strontium compounds, respectively (the su- 
perscript C denotes the values obtained from the analysis 
of the specific heat data). 

To interpret the specific heat data, we compare Cmag 
with the simulated curves for the representative spin- 
1/2 models in both one dimension (uniform chain) and 
two dimensions (non-frustrated square lattice, triangular 
lattice), see Fig. [SJ The maximum value of the mag- 
netic specific heat (C^^g ) and the shape of the max- 
imum are characteristic of the magnitude of quantum 
fluctuations. The enhancement of quantum fluctuations 
suppresses correlated spin excitations, therefore C^^^ 
is reduced, and the maximum gets broader. One may 
follow this effect in Fig. [S] The non-frustrated square 
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FIG. 4; (Color online) Temperature dependence of the specific heat measured at zero field for Ca(VO)2(P04)2 (left panel) and 
Sr(VO)2(P04)2 (right panel). The open circles are the raw data, the solid lines show the phonon contribution Cphon as found 
from the fit to Eq. ([2]), and the dashed lines denote the magnetic contribution Cmag- The insets show the fits with Eq. ((2} in 
the wide temperature range. 



lattice shows a high C^^^ ~ OAGRMi^ In the ID 
case, C^^^ is decreased to 0.35i? due to stronger quan- 
tum fluctuations.— 1^ The triangular lattice (a frustrated 
2D system) shows an even lower (C^fag' — 0.22R) and 
broader maximum as compared to the uniform chaiui^i 

According to the results of band structure calculations 
(Sec. El), the spin systems of M(VO)2(P04)2 are 3D. The 
high dimensionality should result in a narrow maximum 
of the magnetic specific heat with a high absolute value at 
the maximum. However, the M(VO)2(P04)2 compounds 
reveal broad maxima with the absolute values of about 
0.25i? and 0.23i? for M = Ca and Sr, respectively. Such 
absolute values are well below that for the uniform chain 
(ID non-frustrated system) and nearly match the value 
for the triangular lattice (2D frustrated system) . This re- 
sult points to the presence of strong quantum fluctuations 
in the M(VO)2(P04)2 compounds. The spin systems are 
3D; therefore the fluctuations should be entirely caused 
by the frustration. The different breadth of the maxima 
for the calcium and strontium compounds (see Fig. E]) 
may also be an indication of the stronger frustration in 
Sr(VO)2(P04)2. 

The strong frustration evidenced by the specific heat 
data is consistent with the considerable reduction in the 
ordering temperature Tjv as compared to T,^ax (^^^ ^^e 
previous subsection). However, in contrast to the geo- 
metrically frustrated systems, Tjv does not vanish com- 
pletely, and at sufficiently low temperatures, long-range 
magnetic ordering is established. To understand the na- 
ture of the ordered state, we studied field dependence of 
the specific heat (Fig. [6]). At large, the transition tem- 
perature gradually decreases with the increase in the field 
and finally gets suppressed below 0.4 K at 8 and 11.5 T 
for Ca(VO)2(P04)2 and Sr(VO)2(P04)2, respectively^ 



The H — T phase diagrams shown in the insets of Fig. [6] 
are typical for antiferromagnets and point to antiferro- 
magnetic ordering in the M(VO)2(P04)2 compounds. 

C. ^'^P NMR 

For both the compounds, the ^^P NMR spectra con- 
sist of a single and narrow spectral line as is expected 
for / = 1/2 nuclei^ii^^ The single spectral line implies 



chain 




12 3 



FIG. 5: (Color online) Magnetic contribution to the specific 
heat of M(VO)2(P04)2 and theoretical curves for the uniform 
chain (Refs.lsol and lsil) . non- frustrated square lattice (Refs.ls^ 
andfla), and triangular lattice (Ref. [s^) models. The reduced 
temperature scale T/Jc with Jc = yJYlii used. 
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FIG. 6: (Color online) Specific heat of Ca(VO)2(P04)2 (up- 
per panel) and Sr(VO)2(P04)4 (bottom panel) measured at 
different applied fields (H). The insets show the respective 
H — T phase diagrams. 

that both Ca(VO)2(P04)2 and Sr(VO)2(P04)2 have a 
unique '^^P site consistent with the structural data. Ad- 
ditionally, narrow line is also a signature of good sample 
quality. Figure [7] shows the representative spectra for 
Ca(VO)2(P04)2 above the transition temperature. The 
line width and the line shift were found to be temperature 
dependent. 

The temperature dependence of the NMR shift K{T) 
(Fig. [8|) behaves similar to the bulk susceptibility x (see 
Fig. [2]). NMR has an important advantage over x for 
the determination of magnetic parameters. At low tem- 
peratures, extrinsic Curie-type paramagnetic contribu- 
tion often affects the bulk susceptibility, while in case 
of NMR this contribution broadens the spectral line but 
does not contribute to the line shift. Therefore, it is 
sometimes more reliable to extract the magnetic param- 
eters from the temperature dependence of the NMR shift 
rather than from the bulk susceptibility. 

The hyperfine Hamiltonian for '^^P can be written in 
the form H = —jhlA^fS, with I and S being dimension- 
less nuclear and electron spins, respectively. The NMR 
shift K is a direct measure of the uniform spin suscepti- 
bility Xspin- Quite generally, their relation is written as 
follows: 

where S is the temperature- independent chemical shift. 
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FIG. 7: (Color online) ^^P NMR spectra for Ca(VO)2(P04)2 
measured at 4 T by the Fourier transform of the spin-echo 
signal above Tjv. The representative spectra show the line 
shift with temperature. All the spectra are normalized to 
unity by dividing the intensity (J) by the maximum intensity 

(Im ax ) . 



which is almost negligible, Na is the Avogadro constant, 
and Ahf is the transferred hyperfine coupling between 
the electrons and the probing nuclei. Therefore to cal- 
culate Ahf, one should use K vs. Xspm plot with T as 
an implicit parameter. Since above 3 K the extrinsic 
paramagnetic contribution is negligible for both the com- 
pounds, we use the bulk susceptibility x instead of the 
spin susceptibility Xspin- The resulting plots are shown in 
Fig.H The fits yield Ahf = (4920±200) and (5125±200) 
Oe/nB for Ca(VO)2(P04)2 and Sr(VO)2(P04)2, respec- 
tively. The linearity of the K vs. Xspin plots confirms 
that by measuring K{T) we can trace Xspin(r) properly. 
The obtained hyperfine couplings for "^^P are of the same 
order as the values previously reported for vanadium 
phosphates^SiiSiiii and indicate a sizeable hybridization 
of P and V orbitals mediated by 2p orbitals of oxygen. 

To extract the magnetic parameters, we fitted K{T) 
curves above 15 K with the Curie- Weiss law x = C/{T + 
^cw)- The resulting values of the effective moment (/x^) 
and the Curie- Weiss temperature {dew) listed in Ta- 
ble [H These values are in good agreement with that 
obtained from the analysis of the bulk susceptibility. 

For the 1/Ti experiment, the frequencies of the cen- 
tral positions of the corresponding spectra at 70 MHz 
have been excited. For a spin-1/2 nucleus, the lon- 
gitudinal magnetization recovery is expected to follow 
a single exponential behavior. In M(VO)2(P04)2, the 
recovery of the nuclear magnetization after an invert- 
ing pulse can indeed be described by a single expo- 
nential, i (M^)gp) = exp (-t/Ti) + C, where 

M{t) is the nuclear magnetization at a time t after an 
inverting pulse, while Ai and C are time-independent 
constants. Temperature dependences of l/Ti are pre- 
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FIG. 8: (Color online) Inverse of ^^P NMR shift (l/K) vs. T 
for Ca(VO)2(P04)2 (upper panel) and Sr(VO)2(P04)2 (bot- 
tom panel). Circles and squares indicate the experimental 
data, while the solid lines show the Curie- Weiss fits. The 
insets present the K vs. T curves. 



sented in Fig. [TOl The spin-lattice relaxation rates are 
temperature-independent above 20 K and rapidly in- 
crease with the reduction in temperature below 20 K. 
The temperature-independent behavior of 1/Ti is typi- 
cal for the paramagnetic regime with fast and random 
fluctuations of electronic spins.— The critical divergence 
in the vicinity of Tn corresponds to the slowing down of 
fluctuating moments and indicates the approach to the 
state with long-range magnetic ordering. 

Furthermore, the 1/Ti data enable to estimate ex- 
change couplings in the system under investigation. Ac- 
cording to Moriya,^-^ the high-temperature limit of 1/Ti 
for a system of localized moments can be expressed as 
follows: 



Y S{S + 1) V27r 



WE 



k |2 
ij\ 



(5) 



where Af- 

''J 



x,y,z) are the components of the 
hyperfine tensor due to the k^^ magnetic atom. The 
Heisenberg exchange frequency we is defined as uj-e = 
Jc{kB/h)y/2zS{S +l)/3, Jc = V'ET^ is the thermo- 
dynamic energy scale of the exchange couplings, and z = 
Zi = 2 (as introduced in Sec. IIVB|) ,^^ Using the relevant 
parameters and the experimental high-temperature re- 
laxation rates 8 and 8.5 ms~^ (see Fig. [10]), we find J^^ ~ 
13.4 and 14.1 K for Ca(VO)2(P04)2 and Sr(VO)2(P04)2, 
respectively. These values are in reasonable agreement 
with the estimates ~ 10 K from the specific heat (see 
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FIG. 9: (Color online) ^'^P NMR shift (K) vs. bulk suscep- 
tibility (x) with temperature as an implicit parameter. We 
employ the x data collected at the applied field of 4 T, since 
this field was used for the NMR measurements. The solid and 
dashed lines show the linear fits with Eq. (|4]). 



Sec. UVBl) . 

In the spectral measurements at low temperatures 
(slightly above T/v), a broad background signal appears 
along with the central peak for both the compounds. At 
T/v, the central peak vanishes, and the broad background 
signal becomes prominent extending over a large field 
range of about 5 T. This effect is illustrated in the inset 
of Fig. \M for Sr(VO)2(P04)2. The central line disap- 
pears at 1.7 K consistent with the ordering temperature 
of 1.8 K as determined by the specific heat measurements 
(see the inset of Fig. [S] and note that we use the T/v 
value at 4 T, since the NMR measurements are carried 
out at this applied field) . The broadening of our spectra 
upon approaching Tm from above is the usual behavior 
expected for the magnetic ordering. Below Tjv, the spec- 
tral intensity melts into a broad background signal which 
gives strong evidence of a large distribution of internal 
static fields in the ordered state. Similar broadening of 
the spectral line has been observed for many other ma- 
terials in the magnetically ordered state (see, e.g., Refs. 
|4iand|44 



V. BAND STRUCTURE 

The M(VO)2(P04)2 compounds reveal quite compli- 
cated spin systems with numerous superexchange paths; 
therefore the presented experimental data do not allow 
to determine the individual exchange couplings unam- 
biguously. To overcome this difficulty, we turn to band 
structure calculations, estimate individual exchange cou- 
plings and construct a microscopic model of the exchange 
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r(K) 

FIG. 10: (Color online) Spin-lattice relaxation rate (1/Ti) vs. 
temperature (T) for the M(VO)2(P04)2 compounds. The in- 
set shows the field sweep spectra of Sr(VO)2(P04)2 measured 
at the transmitter frequency of 70 MHz and temperatures 
close to Tat. 
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FIG. 11: (Color online) LDA density of states for 
Ca(VO)2(P04)2. The Fermi level is at zero energy. The 
calcium contribution is negligible in the whole energy range, 
therefore it is not presented in the figure. The inset zooms 
the image near the Fermi level, only V 3dxy and oxygen con- 
tributions are shown. 



interactions. As we will show below, the computational 
analysis of the M(VO)2(P04)2 compounds is also rather 
difficult, and one can hardly expect the reliable quan- 
titative estimates of all the exchange couplings in the 
systems under investigation. Nevertheless, we succeed to 
establish a reasonable spin model and provide a plausible 
explanation of the frustration (Sec. IVip . 

Experimental data show the similarity of the magnetic 
properties of the M(VO)2(P04)2 compounds with M = 
Ca and Sr. We calculated band structures for both the 
compounds, analyzed exchange couplings and did not 
find any considerable differences between the two sys- 
tems. Therefore, in the following we will discuss the cal- 
cium compound only. The whole discussion is applicable 
to Sr(VO)2(P04)2 as well. 



A. LDA and tight-binding model 

The LDA density of states plot for Ca(VO)2(P04)2 
is shown in Fig. [TT] Valence bands below —3 eV are 
mainly formed by oxygen orbitals, while the states near 
the Fermi level have predominantly vanadium character 
with an admixture of oxygen. The phosphorous and cal- 
cium contributions to these states are tiny and hardly 
visible in the figure. Note that the energy spectrum is 
gapless in evident contradiction with the green color of 
Ca(VO)2(P04)2. This is a typical failure of LDA due 
to an underestimate of strong electron-electron correla- 
tions in the V 3d shell. Local spin density approximation 
(LSDA)-|-[/ calculations readily reproduce the insulating 
spectrum with an energy gap of 2.5 — 2.8 eV. The latter 
values are in reasonable agreement with the sample color. 



TABLE II: The hopping parameters ti of the tight-binding 
model and the resulting antiferromagnetic couplings jf^^^ for 
Ca(VO)2(P04)2. Only the interactions involving PO4 tetra- 
hedra are listed. The effective on-site Coulomb repulsion po- 
tential UcB ~ 4.5 eV. 





tl 


t2 


ts 


t4 


ts 


t (meV) 


-12 


-33 


36 


-3 


-3 




Ji 


J2 


J-s 


Ji 


Js 


jAPM (j^) 


1.5 


11.3 


13.4 


0.1 


0.1 



To study exchange interactions, we construct an effec- 
tive model that includes the relevant, half- filled bands 
only. Thus, we focus on the bands that are close to 
the Fermi level and have predominantly vanadium char- 
acter. According to the discussion in Sec. [TTl distorted 
octahedral coordination of vanadium gives rise to a non- 
degenerate dxy ground state. Indeed, we find four bands 
formed by V dxy orbitals in the energy range between 
-0.2 and 0.1 eV (see the inset of Fig. [lU and Fig. 
These bands correspond to four vanadium atoms in the 
primitive cell of Ca(VO)2(P04)2 and are used for a tight- 
binding fit of the transfer integrals (hoppings) relevant 
for the AFM exchange interactions. The resulting trans- 
fer integrals (t) are introduced to the extended Hub- 
bard model, and the correlation effects are taken into 
account explicitly via an effective on-site repulsion poten- 
tial Ucff*^ In our case t <C C/cff , hence the Hubbard model 
at half-filling can be reduced to a Heisenberg model for 
the low-lying (i.e., spin) excitations. Thus, we are able to 
estimate AFM contributions to the exchange couplings as 
jAFM ^ 4t2/{7^^^ Similar to Ref.lM we use Ues = 4.5 eV 
- a representative value for vanadium oxides. 
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An advantage of the TB approach is the possibihty to 
estimate all the exchange couplings in the system under 
investigation. Yet the TB fit is sometimes non-unique, 
and this is the case for Ca(VO)2(P04)2 due to the pres- 
ence of numerous NN and NNN hoppings. To get an 
unambiguous solution, we assume that the leading in- 
teractions run via PO4 tetrahedra (i.e., correspond to 
ti — introduced in Sec.[TT|- This assumption looks rea- 
sonable, since magnetic interactions in transition metal 
phosphates are usually mediated by PO4 groups (see, 
e.g., Refs. [2^, m, m, and l46l). Moreover, phosphorous 
gives larger contribution to the states near the Fermi 
level as compared to calcium. Thus, the V-O-P-O-V 
superexchange paths should be favorable. 

The resulting TB fit is in perfect agreement with the 
LDA band structure (see Fig. Table [H hsts the 

hoppings involving PO4 tetrahedra and the respective 
exchange integrals 4i The strongest AFM coupling runs 
between non-parallel structural chains (J3), the other 
strong coupling occurs between parallel chains (J2), and 
all the other AFM couplings are weaker at least by a 
factor of eight. The magnitude of J' s is consistent with 
the experimental data, although the precise values are 
somewhat overestimated. For example, considering J2 
and J3 only, we find Jc = 17.7 K that slightly exceeds 
the J^^ « 13 K estimate from NMR and is well above 
the estimate « 10 K from the specific heat. 



B. LSDA+[/, results 

The spin system formed by J2 and J3 is 3D and non- 
frustrated. To get a frustrated scenario, one has to con- 
sider additionally FM interactions using LSDA-|-t/ calcu- 
lations. The regular approach deals with the construction 
of a supercell enabling different patterns of spin ordering, 
the mapping of the resulting energies onto a Heisenberg 
model, and the estimate of individual exchange integrals 
(see, e.g., Ref. 25). However, this approach seems to be 
inappropriate for Ca(VO)2(P04)2, as the magnetic inter- 
actions are very weak, and the change in the total energy 
due to the variation of spin ordering is of the order of 10 
K (i.e., 10~^ Hartree or ~ 10"^ of the total energy). 
Thus, we have to use the smallest unit cell in order to 
achieve the highest possible accuracy, and only sums of 
the exchange couplings can be estimated. Nevertheless, 
one may expect that the TB results on J^^^ will enable 
to resolve FM interactions (at least, qualitatively). 

The results of the LSDA-f [/ calculations are listed in 
Table [ml The most stable spin configuration is formed 
by ferromagnetic structural chains that are coupled both 
ferromagnetically (for parallel chains) and antiferromag- 
netically (for non-parallel chains). We find the same 
ground state for [/^ = 4 — 6 eV, although the exchange 
couplings strongly depend on the Ud value. This problem 
will be addressed below (Sec. IV C|) after a comparison of 
the LSDA-f [/ results with the experimental data and the 
TB estimates. 



TABLE III: LSDA+U estimates for the exchange integrals in 
Ca(VO)2(P04)2 and the resulting Curie-Weiss temperatures 

^ calc 

Ud (eV) Ji + J2 (K) J3 (K) J4 + J5 (K) ^cw (K) 
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2.9 


4.4 


6.0 


6.7 


5 


-14.0 


4.1 


9.2 


-0.4 


6 


-40.5 


3.7 


14.0 


-11.4 



As we have mentioned in Sec. IIIIl Ud is an adjustable 
parameter that depends on the basis set, crystal struc- 
ture of the compound under investigation, and the ob- 
jective parameters of the calculation. Usually, one has 
to fit Ud to some observable quantity (energy gap, mag- 
netic moments, etc.). It is preferable to use the quan- 
tity that is related to the objective parameters, since 
different objective parameters may require the applica- 
tion of different Ud values. Below, we will employ the 
Curie-Weiss temperature ^cw for this purpose. The 
Curie- Weiss temperature is the second-order term in the 
high-temperature series expansion of the magnetic sus- 
ceptibility. According to Ref. H^, for spin- 1/2 system 

OCW = ^/'iJ2^Z^J^ = '^/'^Y.^J^■ ThuS, 6c^ is dirCCtly 

related to the exchange couplings and provides reason- 
able justification for the choice of the Ud value. 

Clearly, the results for {7d = 5 eV are in good agree- 
ment with the experimental estimate ^cw — —0.3 K, 
while for [/^ = 4 and 6 eV the compliance of the 
experimental and computational results is quite poor 
(see Table IIIip . Thus, the optimal value of Ud for 
Ca(VO)2(P04)2 is 5 eV.^^ Yet in contrast to other spin- 
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FIG. 12: (Color online) LDA band structure of 
Ca(VO)2(P04)2 near the Fermi level (thin light lines) 
and the fit with the tight-binding model (thick darker lines). 
The Fermi level is at zero energy. The notation of k points is 
as follows: r(0,0,0), X(0.5,0,0), Zi(0.5, 0.5, 0), F(0, 0.5,0), 
and ^(0,0,0.5) (the coordinates are given along k^, ky, and 
kz in units of the respective reciprocal lattice parameters 
An/a, An/b, and An/c). The points Z and Z\ are equivalent 
due to the face-centered symmetry of the unit cell. 
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1/2 systems r''^'^ the slight variation oi Ud results in a 
huge change in the exchange couphngs (see Table HIT)) . 
This effect will be further discussed in Sec. IV CI In the 
rest of this subsection, we will focus on the exchange 
couplings calculated with Ud ~ 5 eV. 

According to Tables |TT] and [ml J1 + J2 < J^™+J^™ 
indicating considerable FM contribution to either Ji 
and/or J2. J3 is smaller than J^^^ and positive, hence 
there is also a FM contribution to J3 , but the overall in- 
teraction remains AFM. Finally, J4+J5 > J^™+J^™, 
i.e., the LSDA+J7 and TB estimates for these couplings 
are inconsistent. One may think that additional inter- 
actions contribute to the LSDA-I-J7 value of J4 + J5, 
but such contributions yield 2 — 3 K only^ and do 
not explain the controversy. Our experience shows that 
the TB approach is highly reliable for estimating ex- 
change couplings in spin-1/2 systemsi^i^i^ In contrast 
to LSDA-I-C/, the TB approach includes a simple ex- 
plicit relation between Ji and UcS, while the underly- 
ing LDA calculation is ab initio and does not include 
any adjustable parameters. Therefore, we have to admit 
the failure of LSDA-I-C/ to provide valid quantitative esti- 
mates of the weak exchange couplings in Ca(VO)2(P04)2 
(see Sec. IV Cl for further discussion). Nevertheless, the 
LSDA-|-[/ results arc helpful to get a qualitative under- 
standing of the system under investigation. 

To get an idea about the FM interactions in 
Ca(V0)2(P04)2, we consider vanadium-vanadium sepa- 
rations for Ji — J3. The shortest separation (3.46 A) cor- 
responds to Ji , while J2 and J3 reveal larger distances of 
6.15 and 4.35 A, respectively. FM interactions are short- 
range; therefore it is natural to suggest the strongest FM 
contribution to be that of Ji. Since J^^^ is close to 
zero, the overall interaction Ji should be FM. Unlike Ji, 
J2 and J3 are AFM, although the FM contributions may 
also be non-negligible and reduce the absolute values of 
the exchange couplings. The presence of the ferromag- 
netic in-chain interaction Ji is consistent with our stud- 
ies of other vanadium phosphates having similar chains 
of corner-sharing octahedra.^ 

Thus, the basic microscopic model for Ca(V0)2(P04)2 
includes three exchange couplings: ferromagnetic Ji and 
antiferromagnetic J2 and J3 interactions. The resulting 
spin system is 3D and frustrated (see Fig.[T3land Sec.lVI|). 
Unfortunately, the computational results do not enable 
quantitative estimates of individual exchange couplings 
(see the above discussion and Sec. IV Cp . Therefore, we 
do not make further numerical comparisons with the ex- 
perimental data. Clearly, the presence of both the FM 
and AFM couplings is qualitatively consistent with the 
strongly reduced ^cw- Moreover, band structure calcu- 
lations suggest rather weak (about 10 K) exchange in- 
teractions in Ca(VO)2(P04)2 consistent with the exper- 
imental energy scale Jc = 10 — 15 K. We are convinced 
that the agreement between the band structure results 
and the experimental data is reasonable, especially tak- 
ing into account the complexity of the spin system and 
the weakness of the exchange couplings. 



C. LSDA+[/: The influence of Ud 

To understand the J vs. Ud trends, one should con- 
sider exchange integrals as composed of the AFM and 
FM contributions: J = J^™ + J™. Within one-band 
Hubbard model, J^^^ - t^/U, while J^^ is indepen- 
dent of U. Therefore, the overall J should be reduced 
with the increase in U. The repulsion potential Ud is a 
parameter of the computational method rather than the 
true on-site Coulomb repulsion, and J^f^ is not neces- 
sarily proportional to 1/Ud- Yet the increase in Ud tends 
to improve the localization of 3c? electrons hence sup- 
pressing AFM interactions. Thus, one would expect that 
antiferromagnetic exchange couplings should be reduced 
with the increase of Ud- In fact, the strongest coupling 
usually shows 1/Ud behavior, while other exchange in- 
tegrals are nearly independent of Ud (see, e.g., Refs. [25I 
and [13). If J^™ > J^^, the overall J is antiferromag- 
netic for any reasonable value of Ud- However, for small 
jAFM jFM ^]^g situation will be different, and the re- 
sulting J may be positive for small Ud and negative for 
large Ud- 

The above considerations naturally explain, why FM 
interactions manifest themselves (hence reducing 0cw) 
at J7d > 5 eV only. However, the FM interactions are 
further enhanced with the increase in Ud, and this effect 
has a different origin. LSDA-t-C/ treats electronic cor- 
relations within mean-field approximation; therefore the 
filled states of vanadium are merely shifted to lower ener- 
gies as the repulsive potential Ud is applied. At lower en- 
ergies, oxygen states dominate; hence vanadium-oxygen 
hybridization is enhanced with the increase in Ud- In- 
deed, our calculations reveal the decrease in the mag- 
netic moment of vanadium and the partial spin polar- 
ization of oxygen atoms as Ud is increased. In case of 
Ca(VO)2(P04)2, this effect is even more pronounced as 
compared to other vanadium compoundsi^i^ Moreover, 
the exchange couplings are very weak, and the overesti- 
mated vanadium-oxygen hybridization results in unreal- 
istic results for Ji (note the Ji + J2 value of —40 K at 
Ud = 6 eV). 

In summary, we should point to two complications that 
arise during LSDA-I-J7 calculations for vanadium oxides 
with very weak magnetic interactions. First, one has to 
apply a relatively large Ud in order to suppress AFM 
couplings and to reveal FM contributions. Second, large 
Ud tends to overestimate vanadium-oxygen hybridization 
and to produce unrealistic enhancement of the exchange 
couplings. There is an optimal Ud value [5 eV in case 
of Ca(V0)2(P04)2] that gives rise to a reasonable so- 
lution. However, even at this value both problems are 
retained. In general, we think that LSDA+J7 estimates 
of weak exchange couplings in highly complex structures 
are not accurate enough to provide quantitatively cor- 
rect results. Therefore, we mainly rely on the TB model 
in our analysis. Nevertheless, LSDA-I-J7 estimates facili- 
tate a qualitative understanding of the coupling scenario 
resulting in a reasonable microscopic description of the 
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magnetic behavior. 



VI. DISCUSSION 

Our experimental results show that the complex vana- 
dium phosphates M(VO)2(P04)2 are strongly frustrated 
spin systems with competing FM and AFM exchange 
couplings. Magnetic susceptibility data reveal a max- 
imum at T^ax — 3 K and nearly vanishing Curie- Weiss 
temperatures ^cw < 1 K. For simple spin systems with a 
single leading exchange coupling, T^^^ and ^cw are usu- 
ally close to each other, while the ^cw *C ^max I'egime 
implies the presence of several different magnetic inter- 
actions in the system under investigation. As we have 
stated above, 6'cw is equal to half of the sum of all the 
exchange integrals. Thus, the low ^cw value points to the 
presence of both FM and AFM interactions that compen- 
sate each other. Further on, the low Neel temperatures 
(T/v ~ Tj^ax/^) indicate that long-range ordering in the 
system is impeded by either low dimensionality and/or 
magnetic frustration. 

Low-temperature specific heat provides the informa- 
tion on the spin excitations and presents one of the sub- 
stantial characteristics of quantum magnets. In systems 
with strong quantum fluctuations, the specific heat max- 
imum is reduced due to the suppression of correlated 
spin excitations. The maximum of the magnetic specific 
heat of the M(VO)2(P04)2 compounds (C™fJ) equals to 
0.25i? only. This value is comparable to that for the spin- 
1/2 triangular lattice, a system with strong geometrical 
frustration, and lies well below C^^g typical for non- 
frustrated low-dimensional spin systems. Thus, strong 
quantum fluctuations are present in M(VO)2(P04)2, and 
these fluctuations are likely caused by the magnetic frus- 
tration. 

To reveal the origin of the frustration in the 
M(VO)2(P04)2 compounds, we have performed band 
structure calculations and estimated the individual ex- 
change couplings. Despite the problems discussed in 
Sec. IV CI we succeeded to construct a reasonable micro- 
scopic model that provides a plausible scenario of the 
frustration. The coupling within the structural chains 
(the chains of corner-sharing VOg octahedra, see Fig. [1]) 
is FM (Ji), while interchain couplings are AFM (J2 and 
J3). The interactions form a quadrangle with one FM 
and three AFM sides (Fig. [T31). The FM in-chain inter- 
action favors the parallel alignment of the spins within 
the chain, while the inter-chain AFM interactions favor 
the antiparallel alignment. Thus, the spin system is frus- 
trated. 

Below, we will use several experimental quantities 
{Oc^NtJc, and Hs) in order to estimate the individual 
exchange couplings of the proposed model. Assum- 
ing the ground state revealed by the LSDA-|-[/ calcu- 
lations (see Sec. IV Bl) . we find Hs ~ 2J^g^B/kB (the 
Heisenberg model is treated in a classical way). Us- 
ing the averaged estimates of Jc = [J^ -t- J^^)I2 and 
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FIG. 13; (Color online) Enlargement of the crystal structure 
of Ca(VO)2(P04)2 (see Fig. [ij that shows the proposed sce- 
nario of the frustration. Antiferromagnetic inter-chain cou- 
plings J2, Jz favor antiparallel orientation of spins for neigh- 
boring atoms in the chain (marked by large circles), while 
the in-chain coupling (Ji) is ferromagnetic and would favor 
parallel alignment of the respective spins. 



^cw = (^cw + ^cw)/2 and assuming negative Ji, one 
arrives at Ji = —9.6 K, J2 = 3.5 K, and J3 = 5.4 K 
for Ca(VO)2(P04)2 and Ji = -8.9 K, J2 = 1.8 K, and 
J3 = 7.7 K for Sr(VO)2(P04)2. These values should be 
considered as rough estimates, since our model includes 
three interactions only, and the results of different meth- 
ods for Jc and ^cw are simply averaged. The estimates 
emphasize the similarity of the calcium and strontium 
compounds, although J2 and J3 values are somewhat dif- 
ferent for M — Ca and Sr. The latter result may be rele- 
vant for the different breadth of the specific heat maxima 
in the two compounds (see Fig. [5]). 

We should emphasize that at least two couplings (Ji 
and J3) in the M(VO)2(P04)2 compounds have similar 
magnitude, and the spin system is essentially 3D rather 
than ID with frustrated inter-chain couplings. The frus- 
tration is controlled by the magnitudes of the compet- 
ing exchange interactions. For example, the decrease 
of the absolute value of Ji will reduce the frustration, 
while for AFM Ji the system is non-frustrated at all. To 
the best of our knowledge, no theoretical results for this 
model are available. Basically, one may expect the pres- 
ence of strongly frustrated regions [one is revealed by the 
M(VO)2(P04)2 compounds] and quantum critical points 
in the respective phase diagram. Further studies of the 
model and the appropriate materials are desirable. 
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Finally, we will focus on the structural aspects of 
the present study. The main structural feature of the 
M(VO)2(P04)2 compounds deals with the chains of 
corner-sharing VOg octahedra. Such chains are typical 
for a wide variety of vanadium oxides. ^"^ Within a very 
straight-forward and naive approach, one can identify the 
chains of vanadium polyhedra as spin chains and arrive 
at ID spin system with weak inter-chain couplings. Ac- 
cording to Sec. HU such an approach is inconsistent with 
the orbital state of vanadium, and the actual spin system 
may have any dimensionality depending on the exchange 
couphngs via PO4 tetrahedra or other side groups. In 
the case of M(VO)2(P04)2, the actual spin system is 
3D. Yet in Sr2VO(V04)2 the spin system is ID (con- 
sistent with the naive expectations), but the spin chains 
are perpendicular to structural ones (in contrast to the 
naive expectations).-^ In general, the structures formed 
by VOg octahedra and non-magnetic tetrahedral groups 
(PO4, V+^04, Si04, etc.) provide a promising way for 
studying different spin systems, and further investiga- 
tions could be interesting. However, we should point 
to the importance of careful thermodynamic measure- 
ments and electronic-structure-based microscopic model- 
ing for proper understanding of the respective systems. 
Thus, the interesting physics of the M(VO)2(P04)2 com- 
pounds is well hidden behind the ID features of the crys- 
tal structure and weak magnetic interactions resulting in 
the paramagnetic behavior above 4 — 5 K. Only the com- 
bination of experimental and computational approaches 
enables us to unravel the relevant microscopic mecha- 



nisms and to achieve an insight of the underlying physics. 

In conclusion, we have shown that the M(VO)2(P04)2 
compounds reveal strongly frustrated 3D spin systems 
with competing FM and AFM interactions. The pres- 
ence of both FM and AFM interactions is indicated by 
the vanishing Curie- Weiss temperature ^cWj while the 
strong frustration is evidenced by the reduced maximum 
of the magnetic specific heat. The thermodynamic en- 
ergy scale of the exchange couplings isJc = 10— 15K 
as shown by the specific heat and NMR data. Band 
structure calculations suggest a plausible scenario of the 
frustration with the FM interaction along the structural 
chains and the AFM interactions between the chains. In 
the resulting three-dimensional spin system, the frustra- 
tion is controlled by the magnitudes of the competing 
exchange couplings. The proposed spin model deserves 
further experimental and theoretical studies. 
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